Ultrafast dynamics of a magnetic antivortex - Micromagnetic simulations 
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The antivortex is a fundamental magnetization structure which is the topological counterpart 
of the well-known magnetic vortex. We study here the ultrafast dynamic behavior of an isolated 
antivortex in a patterned Permalloy thin-film element. Using micromagnetic simulations we predict 
that the antivortex response to an ultrashort external field pulse is characterized by the production 
of a new antivortex as well as of a temporary vortex, followed by an annihilation process. These 
processes are complementary to the recently reported response of a vortex and, like for the vortex, 
lead to the reversal of the orientation of the antivortex core region. In addition to its fundamental 
interest, this dynamic magnetization process could be used for the generation and propagation of 
spin waves for novel logical circuits. 



Extended ferromagnetic films often display complex 
magnetization patterns with a rich variety of features 
PP. This complexity can be reduced by decomposing 
the magnetic patterns into a few elementary magneti- 
zation structures, such as domains, domain walls or vor- 
tices. The dynamic properties of such fundamental struc- 
tures have been investigated thoroughly over the last 
years by isolating them in patterned elements In 
particular, the magnetic vortex has attracted much in- 
terest H E]- An equally fundamental, 
yet much less studied magnetic structure is the antivor- 
tex, the topological counterpart of the vortex. In the 
complex structures occurring in extended soft-magnetic 
films, antivortices can be found almost as frequently as 
ordinary vortices: they occur in cross-tie domain walls, 
where they are enclosed by two adjacent vortex struc- 
tures [10] . While the in- plane magnetization distribution 
of an antivortex is very different from that of a vortex 
(see Fig. [TJa), it contains, like the vortex, a tiny core [5] 
at its center in which the magnetization points perpen- 
dicular to the plane. Moreover, the two structures are re- 
lated by underlying topological properties: In both cases, 
the local magnetization rotates by 360° on a closed loop 
around the core. The structures differ by their opposite 
sense of rotation along such a loop, which is quantified 
by the winding number w (w = — 1 for the antivortex, 
w = +1 for the vortex) [TT| [12] . The winding number has 
been predicted to have a direct impact on the magneti- 
zation dynamics [111 113] - However, not much is known 
to date about the dynamic properties of antivortices, 
even though several other fundamental magnetization 
structures have been analyzed in patterned thin-film ele- 
ments [2 J . Studies with high spatial and temporal resolu- 
tion have demonstrated that domains, domain walls and 
vortices exhibit different excitation spectra [5] Q3] . Such 
investigations have mostly focussed on small perturba- 
tions and on reversible changes in these structures pro- 
duced by an external field [15j . However, in a recent 
study on cross-tie walls, Neudert et al. |T7] have reported 
the creation of new cross-ties, i.e., of vortex-antivortex 



pairs, in response to high-frequency magnetic fields. Also 
the field-pulse induced reversal of vortex cores has been 
found to be characterized by the creation of a vortex- 
antivortex pair |?1 U |j. During this core switch, a 
transient cross-tie structure dissolves within a few tens 
of ps via an annihilation process |T2j [18]. Using fully 
three-dimensional micromagnetic simulations we study 
here the ultrafast dynamic response of a single magnetic 
antivortex in response to external field pulses a few tens 
of ps long. We find that a rapid series of pair creation 
and annihilation processes occurs, which is complemen- 
tary to the one observed in a vortex [8] and which leads 
to the reversal of the antivortex core. This represents a 
new fundamental process in magnetism on the nanoscale. 
The spin- wave bursts resulting from the annihilation pro- 
cess make the antivortex structure a possible source for 
spin- waves to be used in novel logical circuits |19j . 

Unlike vortices, it is rather difficult to isolate an an- 
tivortex, i.e., to prepare a nanostructure that contains 
only a single antivortex. This has been achieved by 
Shigeto et al. using thin-film elements consisting of four 
connected rings |20] , where an antivortex was observed 
at the ring junction. Here we use a sample with a 
shape corresponding to the connecting part of four ad- 
jacent rings to stabilize an isolated antivortex (shown in 
Fig. [IJ)). The sample is constructed from four circular 
segments with 200 nm radius. The four short stretch- 
ing strips are 40 nm wide, and the thickness is 20 nm. 
The in-plane shape anisotropy of this particular geom- 
etry can sustain an antivortex due to the tendency of 
the magnetization to align with the sample boundaries. 
The simulations have been performed with the micro- 
magnetic finite-element code that we have used previ- 
ously, e.g., in Ref. |5]. The antivortex structure was ob- 
tained by selecting a suitable initial magnetic configura- 
tion and letting the structure relax [35J . The material pa- 
rameters used in the simulations are those of Permalloy: 
A = 13pJ/m (exchange constant), /j,qM s = LOT (with 
M s the saturation magnetization) , K=0 (anisotropy con- 
stant) and K s — 10 _4 J/m 2 (surface anisotropy) [2"T] . 
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FIG. 1: (color) a: Schematic of an antivortex (left) compared with a vortex (right). The arrows represent the in-plane 
magnetization, b: Modeled concave-shaped sample. The magnetization direction in the branches is schematically represented 
by the grey arrows. The brown area at the center of the structure represents the perpendicular component of the magnetization 
m z , showing the antivortex core, c-j: Distortion of the antivortex structure induced by an in-plane field pulse. The frames are 
snapshots at different times relative to the applied pulse's maximum. The arrows in the top row (c-j) represent the in-plane 
magnetization, while the colors represent the a;-component of the magnetization (m x ). In the bottom row (g-j) the green and 
red "ribbons" are the m x — and m y = isosurfaces, respectively. The colors represent m z . 



About 200,000 tetrahedral elements were used for the dis- 
cretization of the sample. This corresponds to a cell size 
of about 4nm 3 . The magnetization dynamics was cal- 
culated using the Landau-Lifshitz-Gilbert equation with 
damping parameter a = 0.01. The isosurface represen- 
tation introduced in Ref. [T^] has been used to highlight 
the location of the antivortex core, that is, the region 
where the magnetization is exactly perpendicular to the 
sample plane (m z = ±1). This situation occurs where 
the m x = and m y = isosurfaces intersect (m is the 
normalized magnetization: m = M/M s ). 

The antivortex structure was perturbed by a short in- 
planc Gaussian-shaped pulse applied in the y direction. 
In the case shown in Fig. [I] the pulse had a maximum 
intensity of 60 mT and a duration of a — 80 ps. Fig. [l]c- 
j shows the magnetization dynamics in response to the 
applied field pulse. In the top row (c-f), the in-plane 
magnetic structure is shown. The same structure is also 
shown in the bottom row ( g-j), but the colors there repre- 
sent the perpendicular component of the magnetization 
(m z ) and the crossing point of the isosurfaces helps dis- 
play the position of the antivortex core. Following the 
application of the field pulse, the antivortex is displaced 
from its equilibrium position and the in-plane magneti- 
zation of the sample is distorted (d). In a region close to 
the core, this distortion leads to the formation of a "dip" 
in which the magnetization rotates out of the plane in the 
direction opposite to the antivortex core (h). The forma- 
tion of such a dip in the case of an antivortex is analogous 
to the recently reported dip formation near a distorted 



vortex core [5J [53] . Its direction is determined by the 
stray field of the antivortex core [5] . Approximately 80 ps 
after the pulse maximum, a new antivortex is emitted 
from the original antivortex (e). However, according to 
the Poincare-Hopf index theorem of vector fields, the to- 
tal winding number of the magnetization in a thin- film el- 
ement is a topological invariant, and the the formation of 
a single antivortex would violate its conservation. A new 
vortex is thus also produced, which is located between the 
two antivortex structures. The newly formed antivortex- 
vortex pair is made clearly visible by the two additional 
intersections of the isosurfaces (i). This strongly inho- 
mogeneous structure is eventually resolved through the 
annihilation of the initial antivortex with the newly cre- 
ated vortex, whose cores have opposite orientation [12] . 
This results in a sudden reduction of the local exchange 
energy density and leaves behind the newly formed an- 
tivortex with oppositely magnetized core (f, j). 

The production of an antivortex-vortex pair from an 
antivortex structure is a previously unreported micro- 
magnetic process. The resulting transient magnetic con- 
figuration could be called an anti cross-tie wall, i.e. a sin- 
gle vortex enclosed by two anti vortices, which has never 
been observed in a stable magnetic structure. It is also 
remarkable that, in spite of the very different magnetic 
in-plane structure of vortices and antivortices, their ul- 
trafast dynamics are analogous, involving the creation of 
new magnetic structures followed by a destruction pro- 
cess. Here, the strength of the isosurface representation 
is obvious: it captures the main features of the processes 
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FIG. 2: (color) Vortex nucleation at the sample edge dur- 
ing the antivortex-initiated pair creation process. Left: For a 
Gaussian pulse of 60 mT and 80 ps duration, no edge nucle- 
ation occurs. Right: For a stronger pulse (80 mT) of same du- 
ration, the m x and m y isosurfaces cross at the sample bound- 
ary, indicating the nucleation of a vortex. 



unfolding in both vortex and antivortex structures |29j . 
demonstrating the complementarity between them. 

In the case shown in Fig. [l] the system eventually re- 
laxes back to an antivortex structure. However, the sta- 
bility of an isolated antivortex is rather weak compared 
to that of a vortex. This is due to the magnetic flux, 
which in an antivortex has a saddle-point configuration 
while a (energetically favorable) closed flux is obtained 
around an ordinary vortex. The magnetic configuration 
in the sample we have used leads in fact to the formation 
of magnetic surface and volume charges, which reduce the 
stability of the antivortex [3U] . This easily introduces the 
possibility of a vortex nucleating at the sample bound- 
aries if stronger external field pulses are applied. Using 
isosurfaces, the vortex nucleation process at our sample's 
boundaries can easily be described and visualized. The 
behavior of the m x — and m y — isosurfaces is shown 
in Fig. [2] for two field pulse strengths. The crossing of 
the isosurfaces at the boundaries indicates that a vortex 
is nucleated [31]. The introduction of a vortex in the 
sample affects the dynamics of the antivortex and thus 
also the final equilibrium structure. The different possi- 
ble final states are shown in Fig. [3] as a function of the 
applied pulse's strength and duration. It indicates that 
while it is rather difficult to isolate an antivortex, it is 
pretty easy to dissolve it. Starting with low and short 
field pulses, there is a set of parameters for which the 
system is slightly excited, causing the antivortex core to 
rotate about its equlibrium position (region A in Fig. |3|. 
The modes excited by such low field-pulses have been 
investigated in Ref. [15]. When pair production and, 
consequently, a core switch is induced, the studied sys- 
tem relaxes into an antivortex configuration only for a 
rather narrow set of field pulse parameters (region B). 
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FIG. 3: (color online) Final magnetic configurations depend- 
ing on the field pulse strength and duration. The duration is 
defined as the width of the Gaussian. The field was varied in 
increments of 10 ps and for every 10 mT. 



For stronger or longer pulses, vortex nucleation occurs 
concurrently to the antivortex core switch, as shown in 
Fig. [2] The vortex can then migrate through the sample, 
leading to the expulsion of the antivortex after a few ns 
(region C). In region D (for pulses above 60 mT), the 
pulse is sufficiently strong to induce the production of 
new pairs originating from both the antivortex and the 
nucleated vortex. The result is the generation of a tran- 
sient state containing multiple cross-ties. A series of an- 
nihilations [T2J H3 HI] occurs within a few ten ns, leading 
to a a final magnetic configuration consisting only of a 
pair of head-to-head and tail-to-tail domain walls. 

The annihilation of a vortex and an antivortex with 
opposite core magnetization has been shown to be con- 
nected with the emission of spin- wave bursts in the GHz 
range [12j[T8]. It has therefore been suggested that the 
core switch mechanism of a vortex could be used to in- 
ject spin waves into a strip acting as a waveguide |25j . 
The spin waves could then be processed in logical cir- 
cuits [TS]. In this context, the antivortex structure could 
equally be used as a spin-wave source. If its stability can 
be increased, an antivortex would probably be more suit- 
able than a vortex for this purpose: When a disk-shaped 
source is used (in the case of the vortex core switch) a 
90°Neel wall is inevitably present between the source and 
the strip along which the waves propagate. Such do- 
main walls have been reported to have unfavorable prop- 
erties concerning spin wave propagation, mostly reflect- 
ing the spin waves [25]. In contrast, as shown in Fig. [4] 
the spin waves produced in the antivortex can propa- 
gate unhindered into a branch, which naturally extends 
the magnetic configuration of the astroid-likc sample. All 
branches can be extended to connect to other astroid-like 
elements. Such a network could constitute an antidot ar- 
ray |27j and provide a regular lattice of antivortices. 
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FIG. 4: (color) Spin wave propagation at 136 ps (left) and 
226 ps (right) after the pulse's maximum. The spin waves 
generated by the core reversal (left) shortly later propagate 
smoothly into the branch (right). The m z = isosurfaces 
shown in green allow to visualize the wave fronts. The blue 
spot represents the switched core (m z ) using the same scale 
as in Fig. [I] 



What does the study of antivortex dynamics lead to? 
First, the complex and ultrafast modifications that an an- 
tivortex can undergo by applying a short field pulse are of 
fundamental interest. A short field pulse can be used to 
switch the core of an antivortex. This occurs through the 
rapid creation of a new antivortex and the annihilation 
of the original one with a short-lived vortex. Second, 
we find that the static topological complementarity of 
vortices and antivortices is equally exhibited in their dy- 
namical behavior. The dynamics of antivortices deserves 
further investigation since it likely plays an equally im- 
portant role as the dynamics of vortices in the case of 
pulse-induced cross-tie domain- wall transformations |17j . 
In this context it will be necessary to determine whether 
the pair-creation reported in Ref. [17] is directly due to 
the dynamics described above, since larger distances of a 
few hundreds of nm are involved in cross-tie walls. Third, 
the annihilation of transient magnetic structures during 
the antivortex core switching process could be used, e.g., 
for generating spin- waves for logical circuits |19j . 
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